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ABSTRACT 

We report the discovery of HAT-P-15b, a transiting extrasolar planet in the 'period valley', a 
relatively sparsely-populated period regime of the known extrasolar planets. The host star, GSC 2883- 
01687, is a G5 dwarf with V =12.16 . It has a mass of l.Ol±O.O4M , radius of 1.08±0.04.R Q , effective 
temperature 5568 ± 90 K, and metallicity [Fe/H] = +0.22 ±0.08. The planetary companion orbits the 
star with a period P = 10.863502 ± 0.000027 d, transit epoch T c = 2454638.56019 ± 0.00048 (BJD), 
and transit duration 0.2285±0.0015d. It has a mass of 1.946±0.066Mj, and radius of 1.072 ±0.043 Rj 
yielding a mean density of 1.96 ± 0.22 gem -3 . At an age of 6.8^'g Gyr, the planet is H/He-dominated 
and theoretical models require about 2% (10M$) worth of heavy elements to reproduce its measured 
radius. With an estimated equilibrium temperature of ~ 820 K during transit, and ~ 1000 K at 
occultation, HAT-P-15b is a potential candidate to study moderately cool planetary atmospheres by 
transmission and occultation spectroscopy. 

Subject headings: planetary systems — stars: individual (HAT-P-15, GSC 2883-01687) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

Ground-based transit surveys have been very success- 
ful at discovering short-period (P < 5 days) transiting 
extrasolar planets (TEPs) since 2006. Because of the 
short duration of transit events (relative to the orbital 
period), candidates are more likely to be identified by 
higher duty cycle observations. The situation is espe- 
cially demanding toward longer orbital periods, where 
the observational constraint couples with the consider- 
ably lower incidence rates of preferred orbital plane in- 
clinations. The only way we can increase the probabil- 
ity of discovering long-period TEPs from the ground is 
to increase the duty cycle by establishing telescope net- 
works spread in geographical longitude. Surveys such 
as WATTS/TrES and HATNet are designed with this 
goal in mind. For example, in t he particular case of 
WATTS/TrES (jOetiker et al.ll20Toh . for a typical observ- 
ing season, the detection probability of a 10 day period 
TEP from a single observing site is 1-5%. This increases 
to 10% by using two sites of the network. It is clear that 
discovering long-period TEPs with ground-based surveys 
involves a large degree of serendipity even in the case of 
multisite observations. 
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Fig. 1. — Period distributions of the extrasolar planets discovered 
by radial velocity (RADV) method and by photometric transits 
(TEP - a small fraction of TEPs was discovered by the radial 
velocity method and then later confirmed as TEPs). No TEPs are 
included in the distribution of the RADV planets. HAT-P-15b is 
situated at a relatively sparsely populated region in both diagrams. 



In large part due to observational selection effects, 
from the more than 70 TEPsH known today, there are 
only 6 planets with orbital periods long er than 6 days. 
Thre e of them (CoRoT -4b, 6b and 9b, see lFridlund et'alj 
[20101: iDeeeet al.l [2010h were discovered by theCoRoT 
mission; two, HP 17156b (Barbie ri et al.l [20 07). and 



HP 80606b (IMoutou efaLl 120091 iFossev et alj 120091 



Garcia-Melendo fc McCulloughl 120091 and IHebrard et all 



& 



2010T ) came from the photometric followup observations 
of already known planets from radial velocity (RV) sur- 
veys. The only long-period planet originally discov- 
ered by a ground-b ased photometric survey is WASP-8b 
(jQueloz et al.ll2010() . with a period of 8.16 days. The pe- 
riod distribution of TEPs, displayed in the lower panel 

10 See |http://exoplanet.e"u7| as of May 10, 2010 
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of Figure [TJ shows the very sparse population of TEPs 
above 5 days (logP > 0.7). Interestingly, there are also 
relatively few radial velocity (RV) planets discovered in 
the period interval of 10-100 days (top panel). The de- 
pletion of planets in this regime can be explained by 
the stellar-mass-dependent lifeti me of the protoplane tary 
disk and its effect on migration (jBurkert fc Idall2007l and 
ICurridl2009h . 

The 10.86 day period planet HAT-P-15b was discov- 
ered by the Hungarian-made Automated Telescope Net- 
work (HATNet: IBakos et all 120041 ). Operational since 
2003, it has covered approximately 11% of the Northern 
sky, searching for TEPs around bright stars (8 < I < 
12.5 mag). HATNet employs six wide field instruments: 
four at the Fred Lawrence Whipple Observatory (FLWO) 
in Arizona, and two on the roof of the Submillimeter Ar- 
ray hangar (SMA) of SAO in Hawaii. Since 2006, HAT- 
Net has announced and published 14 TEPs. This work 
describes the fifteenth such discovery, the first TEP with 
a period above 10 days discovered by a ground-based 
photometric survey. 

2. OBSERVATIONS 
2.1. Photometric detection 

The transits of HAT-P-15b were detected with the 
HAT-6 telescope in Arizona and the HAT-9 telescope in 
Hawaii. The region around GSC 2883-01687, a field in- 
ternally labeled as G170, was observed on a nightly basis 
between 2005 September and 2006 February, whenever 
weather conditions permitted. We gathered 8538 images 
with an exposure time of 5 minutes and with readout 
time and other operations yielding a 5.5 minute cadence. 
Each image contained approximately 55,000 stars down 
to / ~ 14 mag. For the brightest stars in the field, we 
achieved a per-image photometric precision of 2.5mmag. 

The calibration of the HATNet frames was performed 
by utilizing standard procedures as described in our 
earlier papers. In brief, star detection, astrometry 
and aperture photometry were performed on the cali- 
brated images by using the Two Micron All Sky Sur- 
vey (2MASS) catalog (ISkrutskie et all I2006D and the 
astrometric routine of iPal & Bakos (2006). The re- 
sulting light curves were decorrelated from systematics 
using t he External Param eter Decorrelation technique 
(EPD, IBakos et all 12010ft and the Trend Filtering Al- 
gorithm (TFA^ IKovacs et all [2005) . The light curves 
were searched for periodic bo x-like signals using t he Box 
Least-Squares method (BLS, IKovacs et al.H2002T) . 

We detected a significant signal in the light 
curve of GSC 2883-01687 (HAT-P-15, also known 
as 2MASS 04245952+3927382; a =04 h 24 m 59.54s , 
5 =+39°27'38.3": J2000; V =12.16 : l"Droege et al.l l2006Tl 
As shown in the upper panel of Figure [2J the detection 
of HAT-P-15 is clear, in spite of its long period. The 
occurrence of subharmonic and harmonic components is 
an artefact of BLS (or any other phase-folding period 
search algorithms with arbitrary signal shape) , and indi- 
cates the significance of the detection. During the total 
time span of 135 days, we had three transit events with 
nearly full, and two more with very limited transit cover- 
ages, comprising altogether some 130 data points within 
the transit. Except for one marginally overlapping event, 
all of them were separately caught by HAT-6 and HAT-9, 



stressing the importance of using a network in the search 
for long-period TEPs. 

We searched for additional components in the signal 
by subtracting the best-fitting trapezoidal signal corre- 
sponding to the peak frequency of the above BLS anal- 
ysis. We did not find any component (transit [BLS] or 
Fourier [DFT]) that appeared significant. The null result 
for the DFT search is shown in the lower panel of Fig- 
ure [2] In agreement with the low chromospheric activity 
suggested by the spectroscopic data (see Section I2.3|) we 
can exclude any coherent Fourier signal above ~ 5 mmag. 




V 
T3 

3 



I 



Tj 



I 



Frequency [d ] 

Fig. 2. — Upper panel: BLS frequency spectrum of the TFA- 
filtered HATNet data of HAT-P-15. Lower panel: DFT spectrum 
of the data prewhitened by the best-fit trapezoidal signal derived 
from the BLS analysis above. Both spectra are normalized to 1.0 
at the highest peak. Insets in the panels show, from left to right: 
prewhitening order, peak frequency, transit depth (BLS) or Fourier 
amplitude (DFT) and signal-to-noise ratio. 
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Fig. 3. — Upper panel: Unbinned light curve of HAT-P-15 in- 
cluding all 8538 instrumental / band 5.5 minute cadence measure- 
ments obtained with the HAT-6 and HAT-9 telescopes of HAT- 
Net, and folded with the period of P =10.8635 days. Lower panel: 
Folded/binned light curve, zoomed to the transit. We used 500 
bins in the full period. Error bars show the lcr errors of the bin av- 
erages. The solid line shows the transit model fit in the simplified 
tmdcl & Agol (2002) approximatio n (wi thout limb darkening, see 
ikos et al. 2010) as derived Section 14.21 

Figure [3] shows the folded HATNet light curve, on 
which is superposed the model fit to be discussed in 
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Section 14.21 The flat out-of-transit part of the light curve 
confirms the null result on the presence of additional 
components, in particular of those with periods commen- 
surable with the orbita l periods (secondary eclip se or el- 
lipsoidal variation, see iSirko fc P aczvhskl 120031 ) . which 
might suggest a blended eclipsing binary. The derived 
transit parameters of the HATNet light curve are also 
in good agreement with the considerably more accurate 
followup data to be presented in Section HOI We ob- 
tained a period of P — 10.8645 ± 0.0046 days, a moment 
of transit center of T c = 2453649.974 ± 0.037 [HJD], a 
relative transit duration (first to last contact divided by 
the period) of q = 0.0188 ± 0.0018 and a transit depth 
of 0.0110 ± 0.0008 mag. The HATNet light curve also 
predicts a short ingress time of (0.063 ± 0.065)Pg days, 
again consistent with the followup data. In spite of the 
relatively large la statistical errors quoted, by using the 
HATNet ephemeris we get a difference of only 0.08 days 
between the transit time predicted from the HATNet 
light curve and that observed from the followup data 
(Section RO] and Table [5]) . 

2.2. Reconnaissance Spectroscopy 

The first step in our follow-up investigation of planet 
host candidates is usually to gather low-S/N, high- 
resolution spectra to obtain a basic characterization of 
the host star and exclude obvious blend configurations. 
As in most of our earl ier works, we used the CfA Dig- 
ital Speedometer (DS; ILathamlll992D . mounted on the 
FLWO 1.5 m telescope. For a descr iption of DS fo r recon - 
naissance spectroscopy we refer to lLatham et al.l (|2009f ). 

We obtained 6 DS spectra of HAT-P-15 in 2006 span- 
ning 62 days. The RV measurements showed an RMS 
residual of 0.33kms _1 . Considering that the overall for- 
mal error of the individual measurements is the same 
as the RMS, the observations suggest that there is no 
detectable RV variation above the kms -1 level. This 
rules out a short period stellar companion. Further- 
more, the spectra were single-lined, showing no evi- 
dence of more than one star in the system. Atmo- 
spheric parameters for the star, including the effec- 
tive temperature T c g+ — 5250 ± 100 K, surface grav- 
ity logg*(CGS) = 3.5 ± 0.25, and projected rotational 
velocity v s mi = 1.5 ± l.Okms" 1 , were d e rived as de- 
scribed by [^rriilNeuhIuiiEEGui^he3 dHH) . The 
effective temperature and surface gravity correspond to 
a late G dwarf. The mean heliocentric RV of HAT-P-15 
is +31.23 ± 0.14kms _1 . These stellar parameters, the 
absence of large velocity dispersion and the lack of any 
stellar contamination indicated that the target warranted 
further follow-up studies. 

2.3. High resolution, high S/N spectroscopy 

To obtain high-resolution and high-S/N spectra for 
characterizing the RV variations and to determine the 
stellar parameters with higher prec ision, we used the 
HIRES instrument ()Vogt et al.lll994f) on the Keck I tele- 
scope located on Mauna Kea, Hawaii. We obtained 24 
exposures with an iodine cell and one iodine-free tem- 
plate. The observations were made on 21 nights during 
a number of observing runs between 2007 August 24 and 
2009 December 28. 

The instrumental setup and reduction method used for 
HIRES are the same as given in our earlier papers and 
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Fig. 4. — Top: RV measurements from Keck for HAT-P-15, 
shown as a function of orbital phase, using our best-fit model (see 
Section |2J. Zero phase is defined at the transit midpoint. The 
center-of-mass velocity has been subtracted. Note that the error 
bars include the stellar jitter of 7.17ms -1 , added in quadrature to 
the formal errors given by the spectrum reduction pipeline. Second 
panel from the top: phased residuals after subtracting the orbital 
fit. The rms of the residuals is 7.73ms . Third panel from the 
top: phased bisector spans (BS, mean value subtracted). The low- 
lying scattered points are attributed to the bias caused by the Moon 
(sec Section |3}. Please note that the three most extreme outliers 
have been omitted, in order to show sufficient amount of detail. 
Bottom: Calibrated S activity index values computed from the 
Keck spectra. Note the different vertical scales of the panels. Open 
circles are for the re-plotted values from the [0,1] phase interval. 



described in more detail by M arcv fc Butlerl ()1992f ) and 
iButler et all {1996). The final RV data and their errors 
are listed in Table [1] The folded data are plotted in 
Figured) with our best fit (see Section HJ superimposed. 
Not plotted is the data point taken on BJD 2454370.9769, 
due to its outlier status (most probably due to the large 
contamination from the Moon - see Section [3]) . 

For each spectrum we measured the S and R' HK chro- 
mospheric activity indices from the emission in the Call 
H and K line cores. The S index is computed follow- 
ing the prescription given bv lDuncan et all (|1991| ). after 
matching each spectrum to a reference spectrum using 
a transformation that includes a wavelength shift and a 
flux scaling that is a polynomial as a function of wave- 
length. The transformation is determined on regions of 
the spectra that are not used in co mputing the indices. 
The R' HK index, as described by iNoyes et al.l (|1984| ). 
uses the S-value and the (B — V) color of the star to 
determine the amount of flux in the H and K line cores 
as a fraction of the total energy output of the star. The 
use of (B — V) color in defining R' HK allows for com- 
parison of chromospheric activity across spectral types. 
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TABLE 1 
Relative radial velocity, bisector and activity indices of HAT-P-15. 



BJD RV a ff RV b BS o-bs S c logR' H 

(2,454,000+) (ms -1 ) (ms" 1 ) (ms" 1 ) (ms" 1 ) 

337.11515 -169.94 

339.13623 

339.14209 -85.37 

344.05984 77.73 

345.13121 -3.91 

370.97691 -234.64 

371.97469 -34.16 

372.10061 -32.06 

372.97042 105.93 

373.11653 118.15 

397.91546 115.53 

399.03020 40.57 

427.98281 159.16 

428.91259 167.01 

455.93552 -147.17 

460.87067 160.04 

548.76912 159.50 

728.03347 -179.29 

778.93314 55.05 

809.86948 151.70 

810.84709 116.50 

1109.13606 -192.39 

1192.94775 -32.84 

1192.95593 -23.64 

1193.91529 -87.16 



2.28 


36.70 


6.25 


0.1660 


-4.941 




25.72 


6.78 


0.1670 


-4.935 


2.75 


-14.98 


9.46 


0.1690 


-4.923 


2.40 


26.77 


6.66 


0.1650 


-4.947 


2.35 


25.23 


7.10 


0.1640 


-4.953 


2.76 


-114.84 


15.74 


0.1750 


-4.890 


2.29 


-4.51 


8.90 


0.1690 


-4.923 


2.10 


-57.00 


12.36 


0.1670 


-4.935 


2.15 


2.03 


8.37 


0.1730 


-4.901 


2.03 


7.31 


8.35 


0.1650 


-4.947 


2.75 


-15.73 


10.65 


0.1660 


-4.941 


3.37 


-102.32 


19.47 


0.1750 


-4.890 


2.68 


8.79 


9.32 


0.1620 


-4.965 


3.46 


-17.89 


12.52 


0.1660 


-4.941 


3.07 


29.07 


6.58 


0.1710 


-4.912 


3.45 


26.32 


6.95 


0.1550 


-5.012 


2.44 


30.05 


6.95 


0.1720 


-4.906 


2.33 


22.53 


7.25 


0.1670 


-4.935 


2.66 


27.48 


7.17 


0.1660 


-4.941 


3.05 


26.90 


6.85 


0.1640 


-4.953 


3.53 


28.02 


7.00 


0.1600 


-4.978 


2.72 


-83.90 


14.61 


0.1590 


-4.985 


2.99 


27.35 


6.82 


0.1510 


-5.041 


2.92 


27.56 


6.57 


0.1540 


-5.019 


2.42 


33.33 


6.45 


0.1670 


-4.935 



Note. — The exposure on BJD 2,454,339.13623 was iodine-free, therefore, no RV 
item is given. On BJD 2,454,370.97691 the scattered moonlight affected the derived 
parameters very strongly, including RV. As a result, we left out this observation from 
the analysis. 

a The fitted zero-point that is on an arbitrary scale (denoted as 7 rc i in Section \A.2) 
has been subtracted from the velocities. 

b The values for cr^y do not include the jit ter. 

c Activity indices calibrated to the scale of Duncan e t al.l l|1991l ). 



The Mt. Wilson survey sampled mostly bright, Sun-like 
stars, therefore R' HK is only valid for the color range 
0.4 < B - V < 1.2. Both S and R' HK h ave been cali- 
brated to the Mt. Wilson scale as given bv lDuncan et al.l 
(fl99l . 

The RMS of the residuals of the model fit to the 
radial velocities described later is 7.73ms -1 , which is 
larger than the expected scatter of ~ 2.8ms" 1 , based 
on the formal errors. Assuming that the difference is 
due to spot or other chromospheric activity of the star 
and that the errors add quadratically, we get a value 
of 7.17ms" 1 for the RMS of the "stellar jitter". This 
is more than an order of magnitude larger than the ex- 
pect ed jitter of O^Slnis - 1 , fo r a star similar to the Sun 
(see Makarov ct al. 2009) E3 Indeed, the parameters of 
HAT-P-15 (see Table [3} are very close to those of the 
Sun, albeit its metallicity is slightly higher and its esti- 
mated age is larger, too. The activity indices also imply 
that HAT-P-15 is a very quiet star. The average and 
standard deviation of S and \ogR' HK are (0.165,0.006) 
and (—4.95,0 .04), respectively , matching very nicely the 
solar values (jHall et al.l 120091 ). Furthermore, the low 
chromospheric activity is also in agreement with the sub- 
millimagnitude upper limit obtained for the photometric 
variability (Section l2.ip . This is what we expect for a 

11 However, we note that the average of the measured v alues of 
the stellar jitter in quiet solar-type stars is about 2ms _1 (Wright 
120051) . Nevertheless, this is still a few times lower than the jitter 
of HAT-P-15. 



quiet solar-type star (jMakarov et al.l 120091) . All in all, 
the extra RV scatter is most probably not related to the 
'traditional' source of stellar jitter, but may come from 
other sources, namely: (i) observational errors due to 
instrumental or environmental effects; (ii) light contam- 
ination from the Moon; (iii) perturbation due to other 
bodies in the system. Although the contribution to the 
error budget from (i) may be larger than estimated, here 
we assume that this is not very likely (it would imply a 
factor 2.5 difference between the formal and the true er- 
rors). Possibility (ii) has been examined by omitting the 
data points that are most probably affected by scattered 
moonlight - see Section [3] for demonstrating the strong 
effect of moonlight on the observed distortion of the bi- 
sector span (BS). This test resulted in a decrease in the 
RMS of up to 30%, but only at the price of omitting a 
large number (up to 11) of RV observations that have as- 
sociated BS values affected by the Moon, but otherwise 
fit the single-planet Keplerian orbit fairly well. All of this 
indicates the possible presence of one or more additional 
planets in the system. However, an extensive search for 
the signature of an additional planet in the RV data did 
not yield a significant detection. 

2.4. Photometric follow-up observations 

To confirm the transit signal and obtain high-precision 
light curves for modeling the system, we conducted 
photometric follow-up o bservations with KeplerCam 
(|Szentgvorgvi et alll2005D at the FLWO 1.2 m telescope. 
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5. — Unbinned instrumental z and i band transit light 
curves, acquired with KeplerCam at the FLWO 1.2 m telescope. 
At the bottom of the figure we show the residuals from the fit. 
Error bars represent the photon and background shot noise, plus 
the readout noise. 



TABLE 2 

Journal of photometric follow-up observations of 

HAT-P-15 



0.15 



Date 



N b , 



Cad. [sec] Exp. [sec] Filter 



2007 Dec. 


18 


608 


60 


45 


Sloan z 


2007 Dec. 


29 


514 


15 


30 


Sloan z 


2008 Jan. 


09 


400 


60 


45 


Sloan 2 


2008 Oct. 


17 


276 


75 


60 


Sloan i 



Our journal of observations is given in Table [2] 

The reduction of the images was performed as follows. 
After bias and flat-field calibration, we derived a sec- 
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Fig. 6. — Folded/binned Sloan 2, i band light curve of HAT-P-15. 
The data have been obtained by KeplerCam at the FLWO 1.2 m 
telescope. The barely visible error bars are ltr errors of the bin 
averages. 

ond order astrometric transformation between the ~ 400 
brightest stars and the 2MASS catalog, as described in 
iPal fc Bakosl (pOCJl ). yielding a residual of < 0.2 pix- 
els. Aperture photometry was then performed on the 
resulting fixed positions, using a range of apertures. The 
instrumental magnitude transformation was done in two 
steps: first, all magnitude values were transformed to the 
photometric reference frame (selected to be the sharpest 
image) , using the individual Poisson noise error estimates 
as weights. Next, the magnitude fit was repeated using 
the mean individual light curve magnitudes as reference 
and the RMS of these light curves as weights. In both of 
these magnitude transformations we excluded from the 
fit the target star itself and the 3<r outliers. 

We performed EPD and TEA to remove trends simul- 
taneously with the ligh t curve m odeling (for more details, 
see Section 2] and Bakos et al.1120101 ). In conducting this 
f analysis we reject 3rr outliers from the best-fit transit 
light curve models. For the 2007 December 18, 2007 
December 29, 2008 January 9, and 2008 October 17 ob- 
servations we rejected 4, 4, 1 and 1 images respectively. 
Of the several apertures used each night, we chose the 
one yielding the smallest fit RMS for the light curve. The 
final light curves are shown in the upper plots of Figure[SJ 
To exhibit the resulting photometric accuracy, when all 
data are used, we plot the folded z and i observations 
in Figure [6) Although we plot light curves obtained in 
two different wavebands, there is only a negligible dif- 
ference between the light curves in these two bands due 
to the wavelength dependence of the stellar limb darken- 
ing. The overall standard deviation of the bin averages 
is ~ 0.2 mmag per bin length of ~ 7 min. 

3. EXCLUDING BLEND SCENARIOS 

Here we examine if the distortion of the spectral line 
correlation profiles (characterized by the corresponding 
bisector spans, BS) show any variation correlated with 
the orbital period, radial velocity (or their model-fit 
residuals) and stellar activity indices. It is well known 
that these tests are sensitive to distortions caused by 
hidden background (or physically related) unresolved 
eclipsing binaries or stellar activity (Qucloz et al.l [200TI : 
iTorres et all 12051 . 

Figure 0] shows that the BS distribution as a func- 
tion of orbital phase is nearly flat, except for several low 
points whose deviations are apparently not correlated ei- 
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TABLE 3 
Photometric follow-up of HAT-P-15 



BJD 


Mag a 


"Mag 


Mag(orig) b 


Filter 


(2,400,000+) 










54453.57510 


0.00189 


0.00122 


10.64060 


z 


54453.57987 


0.00257 


0.00121 


10.64190 


z 


54453.58575 


-0.00044 


0.00098 


10.63850 


z 


54453.58644 


-0.00011 


0.00098 


10.63880 


z 


54453.58713 


0.00108 


0.00098 


10.64010 


z 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a The out-of-transit level has been subtracted. These magnitudes 
have resulted from the EPD and TFA procedures, carried out si- 
multaneously with the transit fit. 

b These are raw magnitude values without the application of the 
EPD and TFA procedures. 

ther with the RV residuals or with the activity index. 
Following ou r earlier work related to the verification of 
HAT-P-12b (jHartman et al.l I2009D we investigated the 
contamination caused by the moonlight. By adopting 
the s ame definition o f the s ky contamination factor SCF 
as in lHartman et all (|2009t ). we found a suggestive cor- 
relation between SCF and BS, but with several outliers. 
We attr ibute this to the fact that the simple formula 
given in lHartman et al.l (|2009f ) only holds when the RV 
difference between the Moon and the star is greater than 
the width of the cross-correlation function (CCF) profile. 
When the RV difference is close to zero, the component 
of the CCF due to scattered moonlight will affect both 
sides of the CCF component due to the star, reducing 
the effect of sky contamination on the BS. We addressed 
this by direct modeling of the measured BS. 
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Fig. 7. — Effect of the Moon on the variation of the bisector 
span BS. The model BS values have been derived by synthesiz- 
ing the simplified correlation spectra of the star and that of the 
Moon. Depending on the Moon phase/intensity, target-Moon sep- 
aration and relative moon velocity we get various degrees of con- 
tamination, showing significant correlation with the observed BS 
variation. The apparent outlier corresponding to the measurement 
taken on BJD 2454339.14209 is shown by an open circle. 

For each spectrum we generated an expected CCF, 
taken to be the sum of two Lorentzian functions, one 
for the target and the other for the scattered moonlight. 
The radial velocity separation of the two Lorentzians, 
and their rela t ive in tensities, were determined as in 
lHartman et al.l (|2009f ). These CCFs were used to com- 
pute BS (model) in the same way as it is done for the 



observed CCFs. The result is shown in Figure [7J With 
the possible exception of the measurement taken on 
BJD 2454339.14209, the real BS seems to correlate well 
with the simulated BS, confirming that the former is af- 
fected by moonlight to varying degrees. The cause of 
the one outlier mentioned above is not known, since dur- 
ing the time of the observation the Moon was below the 
Horizon. 

It is also seen that the expected moon-contaminated 
BS variations are a factor of a few lower than the ob- 
served BS variations. Presumably this is due to under- 
estimation of the height of the sky contamination peak 
and/or the use of an incorrect Lorentzian width. Clouds 
will generally increase the sky background beyond the 
estimations, and airglow (which is completely ignored 
here) will also add emission lines to the spectrum at 
RV= 0.0ms" 1 . In the BS analysis we try to clip ob- 
vious airglow lines when cleaning cosmic rays, but this is 
by no means perfect. 

Although visual inspection of Figure 2] suggests that 
there are no correlations between RV, BS and S, we ex- 
amined this problem more rigorously by computing false 
alarm probabilities (FAPs) for various combinations, in- 
cluding data selection based on moonlight contamina- 
tion. Most FAPs are above 30-50% with associated cor- 
relation coefficients CC less than 0.2 (in absolute value). 
The least uncorrelated sets are the RV residuals and the 
BS values, obtained by the omission of the data points 
affected by the Moon (see the 11 points in Figure [7] 
with BS(model) < —5ms -1 ). For the remaining 13 BS, 
RV(residual) data points we obtained a FAP of 3% and 
CC of +0.6. Considering the low number of data points 
in this test, and that the correlation heavily relies on a 
single point, these values are still regarded as characteris- 
tics of a random correlation. It is worth pointing out that 
the relative dispersion of the BS values is fairly small if 
we omit the BS values contaminated by the Moon. The 
ratio of the standard deviation of the RV data (without 
subtracting the orbital motion) and that of the BS points 
is only 0.03. We conclude that all evidence strongly sup- 
port the conclusion that the source of the observed RV 
variation is the orbital motion of a planetary companion 
around HAT-P-15. 

4. ANALYSIS 

For the determination of the physical parameters of 
the planet, first we have to derive those of the host star. 
Since the stellar density is tied to the transit and or- 
bital parameters, we use an iterative process in which we 
utilize several observational and theoretical model com- 
ponents, such as: the high-resolution template spectrum 
obtained with Keck/HIRES; stellar evolution and atmo- 
sphere models; global/simultaneous analysis of the light 
curve and RV data to derive orbital and relative planet 
parameters. We end up with absolute stellar and planet 
parameters and corresponding errors estimated from the 
Monte Carlo simulations performed within the iteration 
process. 

4.1. Properties of the parent star 

We derived the stellar atmosphere parameters from 
the Keck/HIRES template spectrum as given in our 
earlier papers (e.g., iTorres et al.l I2010D and described 
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in more detail i n IValenti fe Piskunovl ([T996) and 
iValenti fe Fischerl (|2005l) . This yielded the following 
initial values and uncertainties^ effective temperature 
Trff* = 5645 ± 80 K, stellar surface gravity log g+(cgs) = 
4.51 ±0.10, metallicity [Fe/H] = 0.27 ± 0.08 dex, and 
projected rotational velocity t>sini = 0.5 ± 0.3kms _1 . 
The differences compared to our previous DS estimates 
(Section 12. 2[) can be attributed in part to fixing the 
metallicity (to that of the Sun) in the DS analysis. 

Following our earlier practice, at this stage we can 
utilize the fact that in the case of a low-mass com- 
panion, the directly measured transit parameters (most 
importantly a/R+, the ratio of the semi- major axis to 
the stellar radius) tightly constrain p*, the mean stel- 
lar density and thereby also logq*, entering in the stel- 
lar atmosphere models (Scager & Mallcn-Ornelas 2003; 
ISozzettiet"aLll2007D . 

Considering the coupling between the orbital and stel- 
lar parameters, in our standard iterative procedure we 
first adopt the values of T c s* , [Fe/H] , and log g+ from 
the initial spectral analysis to fix the q uadratic limb- 
darkening coefficients as given by IClaretl (|2004f ) (these 
coefficients are needed to model the light curves in the 
observed i and z bands). The light curve modeling yields 
a/R+ (see Section |4~2")) . We then use a/R+, T c q+ and 
[Fe/H] to estimate M* and R* by c o mpari son with the 
stellar evolution models of lYi et all (|2001[ ). To obtain 
proper error estimates, the procedure was repeated a 
large number of times with T c fl> and [Fe/H] , drawn from 
Gaussian distributions, whereas a/ R± from the one de- 
rived in the global analysis of Section 14.21 Although the 
resulting surface gravity of log g* = 4.38 ± 0.03 is only 
slightly different from that derived in the initial spectral 
analysis, to see the effect of this change we performed 
a second iteration by holding logg* fixed at the value 
above in the new spectral analysis, and adjusting only 
T eff *, [Fe/H], and vsini. This gave T cSir = 5568 ± 90 K, 
[Fe/H] = +0.22 ±0.08, and vsini = 2.0 ± O-Skms" 1 , 
which we adopt as final atmospheric values for the star. 
The additional stellar parameters have been derived from 
the evolution models (see Table [4]) . 

In Figure [8] we plot the model isochrones of lYi et al.l 
(2001) for [Fe/H]=±0.22 with the final choice of T eff * and 
a/R* marked, and encircled by the la and 2<r confidence 
ellipsoids. 

In supplementing the above stellar parameters, we note 
that the average values of the stellar activity indices 
(see Table [J) also allow us to get independent estimates 
on the stellar age and rotation period. Since HAT-P- 
15 is very similar to the Sun, it is expected that both 
of these parameters ar e very close to those of the Sun . 
Indeed, from Fig. 6 of iMamaiek fc Hillenbrand! (J2008D . 
with \ogR' HK = —4.95 we get an age interval of 4.0- 
8.0 Gyr, in good agreement with our estimate based on 
stellar evolution isochrones (see Table HJ). Similarly, 
the expected rotation rate of ~ 30 days (based on the 
log R' HK calibration of the Rossby number and an as- 
su med overall convec t ive turnover time of ~ 15 day s - 
see iNoves et~aT1 11984 IMamaiek fc Hillenbrand! 120081 ) is 
in good agreement with the value of ~ 27 days derived 

12 Please note that the errors are formal and that we have con- 
servatively increased the error of T e g and [Fe/H] by a factor of 
two. 



from our stellar parameters and spectroscopic v sin i of 
2.0±0.5kms~ 1 . We note, however, that there is no sign 
of significant power of any periodicity in the frequency 
spectra of S and log R' HK ■ 
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Fig. 8.— Stellar isochrones of lYi et al.l l20(jll) in the T eff — 
(a/Ri,) plane for metallicity [Fe/H]=+0.22 and ages of 5.0, 5.5, 6.0, 
6.5, 7.0, 7.5, 8.0, 8.5 and 9.0 Gyr, from left to right. The initial 
value of T e ff * and a/ R* is marked by an open circle, whereas the 
final choice is shown as a solid dot, encircled by the 1<t and 2cr 
confidence ellipsoids. 

As discussed below, HAT-P-15 is projected against the 
edge of the Per OB2 association. Therefore, our usual 
sanity check of confronting the color predictions obtained 
from the evolution models with those measured by stan- 
dard photometry, does not work. Instead, we may get an 
estimate on the reddening and compute the true (dered- 
dened) distance by assuming that the stellar parameters 
(most importantly T c ff) are accurately estimated by our 
spectral and global analyses. The reddening is estimated 
with the help of two, partially overlapping methods. 

In the first method we compare the color index ob- 
tained from the evolution models wit h the one ob - 
tained from the 2MASS survey (jSkrutskie et al.l [20061. 
Here we use the ESO color system and transform the 
J - K inde?0 fr om the 2MASS c olors to the ESO col- 
ors as given by iCarpenteil ()200lD . This yields (J — 
if) (ESO) =0.587 ± 0.033 for the observed color index. 
By using the isochrones value of (Mj - M K ) =0.44 ±0.02 
and the re l ation of E{B - V) = E(J - i<f s )/0.505 from 
IBilir et all (I2008I) . we get E(B - V) = 0.291 ± 0.078. 

A more direct estimate of the reddening can be de- 
rived without the need of a color transformation by us- 
ing the curr ent near infrared color-temp erature trans- 
formation of iGonzalez fc Bonifacio (2009) based on the 
infrared flux method and 2MASS colors. We derive the 
dereddened color index by requiring their formula to give 
the same temperature as the one obtained from our SME 
analysis. At the abundance of ±0.22 ± 0.08 we get 
(J - K s ) = 0.395 ± 0.023, which yields E(B - V) = 
0.313 ± 0.075, in a nice agreement with the one derived 
from the ESO-transformed 2MASS colors. Taking the 
arithmetic mean of the two values, we assign an overall 
reddening value of E(B-V) =0.30 ±0.08 to the system. 

With the above reddening, the apparent magnitude of 
K s =9.641 ± 0.019 and the absolute K s mag of M K = 

13 Unfortunately, direct estimation of E(B — V) is not possible, 
due to the lack of accurate measurements in the B, V filters for 
HAT-P-15. 
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3.13±0.08, derived from the global analysis, we can com- 
pute the true (dereddened) distance of the system. Since 
the interstellar absorption is A(K S ) = 0.382E(B - V) = 
0.115 ±0.029, the true (dereddened) distance modulus is 
6.396 ± 0.087, corresponding to 190 ± 8pc (as usual, the 
uncertainty excludes possible, difficult to quantify, sys- 
tematics which may be present in the model isochrones, 
the color/temperature transformations, or the assumed 
interstellar absorption law). 

Here we examine the relation of HAT-P-15 to the Per 
OB2 association. The average distance, radial veloc- 
ity and proper motion of the association are, respec- 
tively, 318 ± 27 PC +24 ± 4kms~ 1 and 15 ± 1 mas yr~ 1 
(|Belikov et all 120021: iSteenbrugge et all I2003D . From 
our own radial velocity measurements (see Section [2~2j) 
we get that the average radial velocity of HAT-P-15 
is +31.23 ± 0.14km s" 1 . For the proper motioro we 
have ~ 19 ± 6masyr ~ 1 from the UCAC3 catalog (see 
iZacharias et al.l [20101) . We see that the kinematical pa- 
rameters may justify a membership, but our target is 
considerably closer than the mean distance of the asso- 
ciation, although, it is possible that the c luster extends 
some what closer than 200 pc (see Fig. 6 of Belikov et al. 
120021 ). Additional evidence that HAT-P-15 might lie 
within this association is that our analysis above pre- 
dicts a reddening of E(B — V) ~ 0.30, which is close to 
the typical value given by I Belikov et al.l ()2002l ) (see their 
Fig. 4). The distance implies that HAT-P-15 is situated 
in the outer region of Per OB2, nearer to our direction. 
This conclusion is supported also by the rather high value 
of the dust absorption of E(B — V) ~ 0.93 given by the 
map of iSchlegel. Finkbeiner fc Davisl (|1998l ). versus the 
relative low value derived above from the analysis of the 
HIRES template spectrum0 Thus, we conclude that 
HAT-P-15 may well lie within the outskirts of the Per 
OB2 association. Nevertheless we can exclude the possi- 
bility that it is a physical member of the association with 
a high probability, based on its large age of 6.8l 1 ; 6 Gyr. 

4.2. Global modeling of the data 

Our global analysis (involving the light curves and ra- 
dial velocity data, together with some feedback between 
the derived planetary and stellar parameters fully fo l- 
lows the methodology described in Ba kos et al.1 ([2010). 
We summarize here only the most im portant step s and 
refer to the paper mentioned and to iPall ([2009). We 
simultaneously fit the HATNet light curve, the z and 
i band follow-up light curves, and the Keck RV obser- 
vations. Our model for the fo llow-up light curves was 
based on lMandel fc Agoll (|2002l ). while we used the limb- 
darkening free "P1P3" analy tic approximation to this 
model (see iBakos et all I2010D to fit the HATNet light 
curve. The limb darkening coe fficients were computed 
from the tables of lClare t (2004) by using the stellar pa- 
rameters derived from the spectral analysis (Section 14. ip . 

14 The association is oriented in such a way that the projec- 
tion of the true m otion results in a relatively small proper motion 
( StccnbrugEC ct al. 2003|). Therefore, because of the large errors, 
we do not make a comparison between the respective components 
of the proper motions. 

15 We note however that the reddening map of 
ISchlege l, Finkbeiner & Davis ( 1998) may yield somewhat overes- 
timated values at l arge reddenings (cf. [Arcc & Goodman] 119991 : 
ISchuster et"alll2003h . 



TABLE 4 
Stellar parameters for HAT-P-15 



Parameter 



Value 



Source 



Spectroscopic parameters 

T cff * (K) 5568 ± 90 SME a 

[Fe/H] (dex) +0.22 ±0.08 SME 

vsini (kms" 1 ) 2.0 ± 0.5 SME 

Umac (kms" 1 ) 3.70 SME 

« mic (kms" 1 ) 0.85 SME 

7rv (kms" 1 ) +31.23 ± 0.14 DS 

Limb darkening and photometric parameters 

a z 0.2403 SME+Claret b 

b z 0.3173 SME+Claret 

Oj 0.3118 SME+Claret b 

bi 0.3113 SME+Claret 

V (mag) 12.159 TASS 

K (mag,ESO) 9.683 ± 0.020 2MASS+Carpenter d 

Derived parameters 

M^ (M Q ) 1.013 ± 0.043 YY+a/i?*+SME c 

Ri (R Q ) 1.080 + 0.039 YY+a/R t +SME 

logg* (cgs) 4.38 + 0.03 YY+a/iJ*+SME 

U (L Q ) 1.00 ± 0.11 YY+a/iJ*+SME 

M v (mag) 4.86 + 0.14 YY+a/iJ*+SME 

M K (mag,ESO) 3.17 + 0.08 YY+a/iJ*+SME 

M. j - M K (mag.ESO) 0.44 + 0.02 YY+a/_R*+SME 

Age (Gyr) 6.8+^ YY+a/iJ*+SME 

E(B-V) (mag) 0.30 + 0.08 YY+a/.R*+SME 

Distance (pc) 190 + 8 YY+a/R*+SME e 

a SME = "Spectrosco py Made Easy" package fo r analysis of 
high-resolution spectra Valenti & Piskunov (1996). These pa- 
rameters depend primarily on SME, with a small dependence on 
the iterative analysis incorporating the isochrone fit and global 
modeling of the data, as described in the text. 
b SME+Claret = Based on the SME a nalysis and ta bles of 
quadratic limb-darkening coefficients from Claret (2004). 
c YY+a/fl t +SME = YY isochrones fYi ct al. 2001), a/Ri 
luminosity indicator, and SME results. 
d The 2MASS-ESO transformation is based on Carpenter 

G Corrected for reddening 



The transit shape was characterized by the relative plan- 
etary radius p — R p /Ri, the square of the impact pa- 
rameter b 2 , and the reciprocal of the half duration of the 
transit Q/R*. We also included a blending parameter 
-BinstEj for the HATNet data to take into consideration 
the much lower resolution of the HATNet images than 
those taken with KeplerCam. 

As far as the RV data are concerned, we adjusted four 
parameters, corresponding to an unperturbed Keplerian 
orbit: the overall zero point shift 7 rc i of the velocity val- 
ues, the semiamplitude K and Lagrangian orbital ele- 
ments (k,h) = e x (cosWjSincLi). The moment of the 
center of the transit T c is fixed by the photometric data 
and determined by the epoch numbers and by two indi- 
vidual events (of course, assuming strict periodicity rl. 
We assigned the transit number N tr = to the first high 
quality follow-up light curve gathered on 2007 December 
18. The adjusted moments of transit centers in the fit 
were the first transit center observed with HATNet at 

16 ^inst is a scaling factor for the HATNet transit depth: 
^HATNet = B inBt (Rp/Ri) . 

17 We prefer using two epochs rather than an epoch and the pe- 
riod, because we would like to derive quantitie s tha t are as we akly 
correlated as possible - sec Bakos ct al. (2007) and Pal (2008). 
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T Cj _74, and the last high quality transit center observed 
with the FLWO 1.2m telescope, T Ct+2 s- Other events 
were then defined by using these two epochs and the cor- 
responding N tr transit numbers. 

Systematics in the HATNet data were treated by 
EPD/TFA-filtering prior to the global analysis. For the 
FLWO data we employed a joint iterative filtering of the 
above type together with the search for the best fitting 
physical model. The five EPD parameters were the hour 
angle, the square of the hour angle, and the stellar pro- 
file parameters (equivalent to FWHM, elongation and 
position angle). The EPD procedure was performed in 
"local" mode with separate EPD coefficients defined for 
each night. The five EPD parameters were augmented 
by the out-of-transit magnitude of the individual events. 
The TEA filtering was performed in "global" mode using 
the same set of stars and TEA coefficients for all nights. 
We chose 20 template stars that had good quality mea- 
surements for all nights and on all frames, implying an 
additional 20 parameters in the fit. The number of fitted 
parameters were the following: 4 x 6 = 24 EPD, 20 TFA, 
3, related to instrumental configuration (the blend fac- 
tor i?i ns t and the out-of-transit magnitude, Mo,HATNet of 
HATNet and the relative RV zero-point 7 ro i) and the 8 
physical model parameters (transit times T c ._74, T c+ 28 
- for transit epoch and period, R p /R±, b 2 , (/R±, K, 
k = ecoscj, h = esino;). The total number of fitted 
parameter is 55, which is well below the number of avail- 
able data points (1786 photometric followup and 23 ra- 
dial velocity data) . For additional technical details see 
iBakos et all (|2010T) and EH (|2009t) . We note only that 
the errors were obtained through Monte-Carlo simula- 
tions as described, e.g., by Ford (2006). 

The final solution for the relevant planetary parame- 
ters are summarized in Table [5r^ 

Also listed is the "RV jitter" , which is a component of 
assumed astrophysical noise intrinsic to the star that we 
add in quadrature to the RV measurement uncertainties 
in order to have % 2 /dof = 1 from the RV data for the 
global fit (see also Section 12.31 for further discussion of 
the jitter). For possible future observation of occulta- 
tion events we computed the eclipse parameters. Since 
the transit happens near apastron, the occultation occurs 
with certainty. 

5. DISCUSSION 

The transiting extrasolar planet (TEP) HAT-P-15b re- 
ported in this paper is among the few with orbital periods 
greater than 8 days 13 We summarize the properties of 
the currently known long-period TEPs in Table [6] 

HAT-P-15b is the longest period TEP discovered by 
ground-based photometric surveys. With a duration of 
5.5 hours, the transit of HAT-P-15b is at the limit in 

18 Some auxiliary parameters (not listed in Table [5j 
are: T c _ 74 = 2453649.98149 ± 0.00217 (BJD), T c , +28 = 
2454758.05871 ± 0.00073 (BJD), 7 rcl = -14.4 ± 2.2ms"' 1 . This 
average velocity for the Keck RVs is merely a zero point value 
and does not correspond to the true center of mass RV of the 
system (for this see the corresponding DS value in Table [3), 
Binstr = 0.95 ± 0.04. 

19 Currently, there are no TEPs between 5.6 and 8 days, and, 
together with HAT-P-15b, there are only four TEPs with periods 
greater than 10 days. The periods are more densely distributed 
below 5 days. 



TABLE 5 
Orbital and planetary parameters of HAT-P-15 b 



Parameter 



Value 



Light curve parameters 

P (days) 10.863502 ± 0.000027 

T c (BJD) a 2454638.56019 ± 0.00048 

Ti4 (days) b 0.2285 ± 0.0015 

T12 = T34 (days) b 0.0237 ± 0.0014 

a/R* 19.16 ±0.62 

CAR* 9.76 ±0.03 

R p /R* 0.1019 ±0.0009 

b 2 H7+ - 047 

6 = acosi/-R* 0.342l°°86 

i(deg) 89.1±^ 

RV parameters 

K (ms" 1 ) 180.6 ±3.5 

fc RV c -0.025 ± 0.005 

h KV c -0.188 ± 0.019 

e 0.190 ±0.019 

U) 262 ± 1° 

o-Rv(fit) (ms" 1 ) 7.73 

RV jitter (ms" 1 ) 7.17 

Secondary eclipse parameters 

T s (BJD) 2454643.816 ± 0.039 

T s 14 (days) 0.1607 ± 0.0066 

Ts',12 (days) 0.0156 ± 0.0009 

Planetary parameters 

M p (M j) 1.946 ± 0.066 

R p (Rj) 1.072 ± 0.043 

C(M P ,R P ) d 0.48 

p p (g cm -3 ) 1.96 ± 0.22 

a (AU) 0.0964 ± 0.0014 

logg p (cgs) 3.62 ± 0.03 

T cq (K) 904 ± 20 

C 0.346 ±0.015 

Fper (10 8 ergs _1 cm -2 ) £ .... 2.25 ± 0.205 

F ap (lOSergs^cm- 2 ) f 1.05 ±0.11 

(F) (lO^rgs- 1 cm" 2 ) f 1.51 ± 0.137 

a All Julian dates in this paper are based on the UTC scale. For 
applications requiring precisions of the order of a second, such 
as the study of transit timing varia tions, the use o f Terrestrial 
Time is recommended (see, e.g. Torres et al. (2010)). 

T14: total transit duration, time between first to last con- 
tacts; T12 = T34: ingress/egress duration, time between first 
and second, or third and fourth contacts. 

c Lagrangian orbital parameters derived from the global mod- 
eling, and primarily determined by the RV data. 
d Correlation coefficient between the planetary mass M p and 
radius R. p . 

The Safronov number is given by © = i(Vesc/V or b) 2 = 
(a/Rp)(M p /M ic ) fsee lHansen fc Barman]|2007r) . 

1 Incoming flux per unit surface area. 

length that can be pursued from a single site ground- 
based observation (the only planet with a longer period 
than HAT-P-15b that can still be followed up from a 
single site is HD 17156b). The possibility of performing 
single-site full followup observations is a valuable prop- 
erty, since systematics are more difficult to handle when 
low signal-to-noise ratio data from different sites are com- 
bined. The orbit is eccentric with a high significance. 
Except for CoRoT-4b (and perhaps for CoRoT-6b), all 
of the 7 long-period planets have this property. Most of 
the planet masses are greater than 2Mj. The host star 
masses are lower than ~ 1.2M Q , in agreement with the 
upper mass limit of ~ 1.4M© of the extrasolar planets 
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TABLE 6 
Long-period TEPs a 



Name 



M p 
[M 3 



Rp 

[Rj] 



J'-icor 



rp C 

[K] 



Period a e T 14 , Star(T cft ) [Fc/H] Age c M* R* V 

[d] [AU] [h] [K] [dex] [Gyr] [M Q ] [R Q ] [mag] 



WASP-8b 


2.25 


1.05 


52 


940 


8.15872 


0.080 


0.31 


3.5 


G6 (5600) 


+0.17 


1.0 


1.03 


0.95 


9.9 


CoRoT-6b 


2.96 


1.17 





1020 


8.88659 


0.086 


0.10 d 


3.8 


F5 (6090) 


-0.20 


3.0 


1.06 


1.03 


13.9 


CoRoT-4b 


0.72 


1.19 





1070 


9.20205 


0.090 


0.00 


1.1 


F0 (6190) 


+0.00 


1.0 


1.10 


1.15 


13.7 
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1.95 


1.07 


10 


910 


10.8635 


0.096 


0.19 


5.5 


G5 (5568) 


+0.22 


6.8 


1.01 


1.08 


12.2 


HD 17156b 


3.21 


1.02 


136 


880 


21.2169 


0.162 


0.68 


3.2 


GO (6079) 


+0.24 


3.1 


1.24 


1.45 


8.2 


CoRoT-9b 


0.84 


1.05 


18 


110 


95.2738 


0.407 


0.11 


8.1 


G3 (5625) 


-0.01 


3.0 


0.99 


0.94 


13.7 


HD 80606b 


4.08 


0.98 


195 


400 


111.436 


0.455 


0.93 


11.9 


G5 (5574) 


+0.43 


7.6 


1.01 


1.01 


8.9 



a All parameters are from http://exoplanet.eu/ except for those of WASP-8b (sec Qucloz et a].,2010) and HD 80606b fsee lMoutou et all 
120091 (T.ff) and lHebrard et al. 2010). 

Based on the models of lFortnev et all (1200 71 ): zero core mass indicates that the planet has larger radius than predicted by the models at 
zero core mass. 
c Computed at a distance given by the semi-major axis. 

Upper limit. 
e Ages have been taken cither from 'exoplanet.eu' or from currently published papers. For CoRoT-6b we adopted an intermediate value 
from the range of 2.5-4.0 Gyr listed bv lFridlund et al.l ( 1201011 . For CoRoT-9b IDeeg et ahl (1201011 give a wide range of 0.2-8 Gyr, and we 
adopted a plausible intcrmediate-to-low age. 



situated in the 'period vall ey' between ~ 0.1 AU and 
- 1 AU (IBowler et al.ll2009l) . 

The relative depletion of planets in the 0.1-1 AU re- 
gion was recognized a few years ago when the number 
of extrasolar planets grew to the level t hat made more 
reliab le statistical investigations possible ([Johnson et al.l 
2007). All host stars with semi- major axes smaller than 
~ 1 AU have masses lower than ~ 1.4M Q . Above 
~ 1 AU there are several stars with M ~ 0.5-2.5Af Q 
(jBowler et al.| [2009). The distribution of the periods (or 
semi-major axes) of all known extrasolar planets is bi- 
modal with a high significance (see upper panel of Fig- 
ure [lj . Although the sample of multiple systems is not 
very extensive, it is interesting to note that the distribu- 
tion bec omes more un i form w hen considering only these 
systems iWright et al.l (|2009f ). Although some aspects 
of these properties can b e explained by engulfm ent of 
planets in evolved giants (jVillaver fc Livid I2009T ). early 
stage of planet formation should play more important 
role. Acco rding to the models co nsidering the forma- 
tion phase (Burkcrt & Ida 2007 and Curric 2009, see also 
iKretke et al. 2009 for the role of magnetorotational insta- 
bility in the protoplanetary disk) the bimodal distribu- 
tion is attributed to the stellar-mass-dependent lifetime 
of the protoplanetary gas disk. For higher mass stars, the 
gas-depletion time scale is shorter due to accretion and 
increased EUV radiation, and therefore planets may stop 
inward migration too early and get stranded on orbits at 
high semi- major axes. For stellar masses < 1.2M Q plan- 
ets do migrate. The rate of migration is determined by 
the ratio of the viscous and depletion time scales Pi 

We investigated the agreement between the derived 
planetary mass and radius and predictions from cur- 
re nt theoretical model s. We interpolated the models 
of lFortnev et a l. (2007) to the observed planet mass of 
1.946 ±0.066 Mj and semi-major axis of equivalent solar 
irradiance (which is the same as the observed semi-major 
axis in this case). The iso-core lines in the log(age)- 
log(i?p) plane are shown in Figure^ Although the errors 
are large, and the derived core mass (especially at the 

20 With a careful filtering of the currently available data we 
found that the valley exist at all host star masses, albeit for stars 
with > 1.2M0 the effect is indeed much stronger. 



low end) depends very sensitively on the stellar/planet 
parameters, it seems likely that HAT-P-15b is a H/He- 
dominated gas giant planet with a core mass most likely 
about lO-M^ , but which could range from zero to per- 
haps 50M e . 

HAT-P-15b does not exhibit the "radius anomaly", 
where - as is frequently the case for hot Jupiters - the 
planetary radius is signific antly larger than pr edicted 
by models such as that of iFortnev et al.l (|2007l) . Fur- 
thermore, we find that for a hypothetical planet with 
the same mass and radius as those of HAT-P-15b, the 
derived core mass is only a weak function of the pe- 
riod (semi-major axis) between ~ 10 and 1000 days (it 
varies between and 10M®). However, below 10 days 
the core mass increases up to 60M ffi (at 1 day), and 
starts especially being sensitive below 4 days. These all 
show that HAT-P-15b resides in the period regime where 
standard planet structure models (i.e., the ones already 
tested on Solar System planets) are satisfactory, without 
resorting to additional effects (e.g., extra heat source, 
tidal heating, enh anced atmospheric opacities, etc., see 
iMiller et al.ll2009l ) required by many of the hot Jupiters. 

In Table |6] we also added the derived core masses for the 
other long-period TEPs. We see that two of the shorter 
period members still exhibit some radius an omalies. For 
CoRo T-4b and CoRoT-6b, tidal heating ([Miller et al.1 
2009) is unlikely, since they both have close to circu- 
lar orbits. On the other hand, considering the possi- 
ble ambiguities in the observed absolute planetary radii, 
it is possible that future followup observation will re- 
sult in downward radius corrections of 3-5% necessary to 
solve the apparent anomaly for these planets. This level 
of radius correction is not uncommon in the followup 
works of increasing accuracy (e.g., the case of HAT-P-lb 
-iWinn et al.ll2007l or that of HP 80606b - iFossevet all 
120091 versus iHebrard et"aT1l2010L and especially WASP- 
8 b, with 10% radius decrease in the very recent analysis 
of lQueloz et al1l27)Tol ). 

It is also useful to compare the overall dependence of 
the planetary radius on the incident stellar flux. Fig- 
ure[T0]shows this dependence for TEPs in the mass range 
of (0.7-5.0)Mj. Although the stellar fluxes for the long- 
period TEPs are still 2-3 orders of magnitude larger than 
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that of Jupiter, it seems that this is already enough to 
cause a substantial decrease in the planetary radius, and 
bring down the otherwise inflated values to 'more nor- 
mal' ones, corresponding to that of the Jupiter. We note 
that this happens in a quite wide mass range. 




3 3.5 

log(age) [Myr] 

Fig. 9. — Position of HAT-P-15b on t he ag e — planet radius 
diagram. The models of Fort nev et al.l 1)20071 ) have been inter- 
polated to the derived mass and solar equivalent semi-major axis 
of H AT-P-15. Open circl es are the tabulated/interpolated values 
from[Fbrtncv et al. (2007), continuous lines have been obtained by 
linear interpolation. Lines are labeled by the core mass (relative 
to the total planet mass). 




10000 100000 1e+06 1e+07 1e+08 1e+09 1e+10 1e+11 
<F> [erg s" 1 cm" 2 ] 

Fig. 10. — Incident stellar flux versus planet radius for TEPs of 
0-7 < Afp/Afj < 5.0. The long-period planets of Tableware shown 
by whiter shade (green), the triangle at the left corner shows the 
position of Jupiter. 

With V=12.2 mag, HAT-P-15 is still bright enough to 
be the subject of ground- and space-based observations. 
At a distance of a = 0.0964 ± 0.0014 AU from its host 
star, HAT-P-15b allows us to investigate the properties 
of a colder planet, not influenced by the complicated hy- 
drodynamical and radiative effects which may apply to 
canonical hot Jupiters (see, e.g., the discussion of the 
core mass insensitivity above). At the same time, the 
period of 10.9 days is still on a relatively short time scale, 
allowing frequent followup observations. The feasibility 
of possible follow-up observations is briefly summarized 
below. 

• Search for other planets. Long-term radial velocity 
(RV) and transit timing monitoring is justified due 
to the signature of excess RV scatter that is un- 
likely to be due to stellar activity (see Section [23)1 . 
Although we have made a very thorough search 
both in the RV and in the photometric data for a 



trace of a second planet, we have not found any- 
thing convincing. Due to the relatively small value 
of 7.17m s^ 1 of the "jitter", the mass of the pos- 
sible hidden companion is estimated in the sub- 
Neptune regime (except if there is a planet in, e.g., 
a 2:1 resonance, hiding in the e ccentric solution 
- see Anglad a-Escude et al.l l2010h . We note that 
here, the eccentricity of the orbit of HAT-P-15b is 
a less forcing argument for searching for a putative 
second planet, because at the semi-major axis of 
HAT-P-15b the circularization time scale due to 
tidal dissipation may be on the order of several 
Gyrs (jMatsumura et al.ll200cf ). 

• Inclination of the orbital and stellar spin axes. 
The measurement of the Rossiter-McLaughlin (R- 
M) effect for HAT-P-15 is somewhat challenging, 
due to its relatively low projected rotation rate of 
v sin i=2.0 ± 0.5 . The expected size of the effect is 
~ 20m s _1 , similar to the one curre ntly measured 
for WASP-8b with the same vsmi (jQueloz et al.l 
I2010D . Three of the seven long-period planets listed 
in Table [S] have yielded R-M data sufficient to com- 
pute the projected inclination of the stellar ro- 
tational and orbita l axes. HD 17156 is aligned 
(iNarita et all 120091) while HD 806 06b is tilted 
Hebrard et al.1120100 . and WASP-8b (IQueloz et al.1 



2010) is apparently strongly tilted, suggesting a ret 
rograde orbit. It is clear that measuring the R-M 
effect of long-period planets in more evolved sys- 
tems such as HAT-P-15 is very important, due to 
the expected relaxation of the dynamics to a final 
equilibrium state. 

• Thermal emission. With a proper choice of the 
infrared waveband and repeated observations, the 
detection of the occultation of HAT-P-15b by 
ground-based instruments may be feasible. As- 
suming blackbody radiation, at a fixed wave- 
length the occultation depth is equ al to <5 OC c = 
(i?p/-Rstar) 2 (T pl /T star ) (IW^ l20T0l). Assuming 
T cq = T p i, at the moment of occultation (which is 
the periastron in the case of HAT-P-15b) we have 
T p i - 1000 K. It follows then that <5 occ - 0.002, 
which is measurable, especially if we consider that 
the occultation lasts for ~ 4 hours; short-enough 
to cover the event in one night and long-enough to 
gather sufficient number of data points (especially 
in the case of repeated observations) to reach the 
level of detection of several a. 

• Atmospheric absorption. Measuring the small 
change in the transit depth due to the varying ab- 
sorption level in different wavelengths (i.e., per- 
forming transmission spectroscopy) is even more 
challenging than detecting the occultation event. 
This is because the change in the transit depth is 
proportional to the scale-height H, which is pro- 
portional to the rati o of the loc al temperature 
to the gravity T/g |Wmn1 [20Tol) . Since HAT- 
P-15b is relatively cold and has a higher grav- 
ity than 'standard' hot Jupiters, it is expected 
that the signal will be about 25% of what is usu- 
ally measured (assuming 'standard' hot Jupiter T 
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and g of 2000 K and 20m s 2 ) . Considering that 
the relative radius variation even for short-period 
targets like HP 189733 is of the order of 0.1% 
(Sing et al.l I2009T ). measuring absorption features 
in colder planets such as HAT-P-15b would indeed 
be very difficult. 

With the continuing work of ground-based wide-field 
photometric surveys, we expect to discover further TEPs 
in the 10-30 day period regime. Since these will be 
mostly bright (V < 13 mag) targets with orbital periods 
still considered to be relatively short, a wide range of 
followup works will be feasible. Therefore, these planets 
will fill the gap between the classical hot Jupiters and 
the long-period, but likely much fainter planetary sys- 
tems (more akin to the Solar System) to be discovered 
by ongoing and future space projects. 
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